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The temporal and spatial components of the heavy-light vector current and the spatial compo- 
nents of the axial current are expressed in terms of lattice-regulated operators suitable for simu- 
lations of B and D mesons. The currents are constructed by matching the appropriate scattering 
amplitudes in continuum QCD and a lattice model to one-loop order in perturbation theory. In 
the lattice theory, the heavy quarks are treated using the nonrelativistic (NRQCD) formulation and 
the light quarks are described by the tadpole-improved clover action. The light quarks are treated 
as massless. Our currents include relativistic and discretization corrections through 0(a a /M, aot s ), 
where M is the heavy-quark mass, a is the lattice spacing, and a s is the QCD coupling. As in our 
previous construction of the temporal component of the heavy- light axial current, mixing between 
several lattice operators is encountered at one-loop order, and 0(aa s ) dimension-four improvement 
terms are identified. 

PACS number(s): 12.38.Gc, 12.39.Hg, 13.20.He, 14.40.Nd 



I. INTRODUCTION 



An important goal of lattice gauge theory is to pro- 
vide estimates of the hadronic matrix elements, such 
as those of the electroweak currents and effective four- 
fermion operators, which are needed for precision tests 
of the standard model. Numerical simulations of quarks 
and gluons in lattice-regulated QCD currently provide 
the only means of calculating such matrix elements from 
first principles taking long-distance QCD dynamics fully 
into account. An important ingredient in these calcu- 
lations is the construction of the lattice current oper- 
ators which match the currents defined in the contin- 
uum to some desired accuracy. In this paper, we ex- 
press the heavy-light vector and axial-vector currents in 
terms of operators defined in a lattice model in which 
the heavy quarks are treated using the nonrelativistic 
(NRQCD) formulation [jlj, the light quarks are described 
by the tadpole-improved clover action Q], and the glu- 
ons are governed by the standard Wilson action. These 
heavy-light currents are important in studies of heavy 
meson leptonic and semilcptonic decays (for recent re- 
views, see Refs. [§-|[). The lattice currents are deter- 
mined by matching scattering amplitudes in the lattice 
model to those in continuum QCD to one-loop order in 
perturbation theory. The matching is carried out through 
0(1/M, a s /M, a a s ) where M is the heavy quark mass, a 
is the lattice spacing, and a s is the QCD coupling. The 
light quarks are treated as massless. 



Each heavy- light current J M , defined in some contin- 
uum renormalization scheme, can be written 

■h = E C i J Sat + 0(a 2 s , a 2 , l/Ml a s a/M ), (1) 



where jjf} &t are operators defined in the effective lattice 
theory, Mq is the heavy quark (bare) mass parameter ap- 
pearing in the lattice NRQCD action, and the Cj 7 coeffi- 
cients are c- numbers which depend only on a s and aMo- 
The goal is to identify the operators lat and to cal- 
culate the matching coefficients Cj (a s ,aMo). The pro- 
cedure for doing this was described in a previous article 
Q in which the temporal component of the axial-vector 
current was studied. For convenience, we reiterate the 
salient steps: (1) select a quark-gluon scattering process 
induced by the heavy-light current of interest and calcu- 
late the one-loop amplitude for this process in continuum 
QCD; (2) expand the amplitude in powers of 1/M; (3) 
identify operators in the lattice theory, usually by inspec- 
tion, which reproduce the terms in this expansion; (4) 
calculate the one-loop mixing matrix of these operators 
in the lattice theory; and (5) adjust the Cj 7 coefficients 
to produce a linear combination of lattice current oper- 
ators whose one-loop scattering amplitude agrees with 
that determined in step 2 to a given order in 1/M and a. 
In this paper, we apply the above procedure to the spa- 
tial components Ak of the axial-vector current and the 
spatial and temporal components of the vector current 
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VJj . We omit much of the calculational details since they 
are similar to those already described in Ref . |J . 

In Sec. ||, we start from a continuum QCD calculation 
and identify the current operators in the effective the- 
ory needed to reproduce the continuum current through 

0(a s /M). In Sec. Ill, we describe the lattice theory and 

— • ( i) 

the one-loop mixing calculation among the J^L- The 

matching procedure is then completed in Sec. [fVL In 

Sec. [v| we discuss 0(a a s ) corrections in the static limit, 

and issues pertaining to terms that behave as a s log(aM) 

are dealt with in Sec. VI. The paper concludes with a 



summary of our results in Sec. VII 



II. OPERATOR IDENTIFICATION 

The heavy-light vector and axial-vector currents are 
given as usual by V^(x) = q(x)j fl h(x) and A^(x) — 
q(x) 757^ h(x), respectively, where q(x) is the light quark 
field, h(x) is the heavy quark field, and 7^ are the stan- 
dard Dirac 7-matrices in Euclidean space-time which sat- 
isfy {7^,7,,} = 2<^„, and 75 = 7071 7273- Euclidean- 
space four-vectors are defined in terms of Minkowski- 
space four- vectors (indicated by an underline) using xq = 



and x 



3 — 



-Xj, for j = 1, 2, 3. For the derivative 
operator, <9 = — id and dj = dj = —& ■ The currents are 

related by Jq = Jx), and J- 7 ' = —iJ_ J . For further details, 
see Refs. §0. 

The first step in identifying the necessary operators in 
the effective lattice theory is to calculate the matrix el- 
ements of the continuum QCD currents for an incoming 
heavy quark of momentum p and an outgoing light quark 
of momentum p' . Using the on-shell mass and wave func- 
tion renormalization scheme in Feynman gauge and ex- 
panding in 1/M (except for Uh(p)), one finds to one-loop 
order in perturbation theory 



(q(p')\J»\h(p)) QCD 
where 



u q (p')W'(p',p)u h (p), (2) 



W 3 u {p',p)=a l rr 
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For the vector current, TY = 7^ and for the axial- vector 
current, = 757^. Uh(p) and u q (p') are the stan- 
dard spinors for the heavy and light quarks, respectively, 
which satisfy the Dirac equation. The light quark mass is 
set equal to zero. Ultraviolet divergences are treated us- 
ing dimensional regularization, and fully anti-commuting 
75 matrices are used. We use a gluon mass A to regulate 
infrared divergences. Note that M is the heavy-quark 
pole mass. 

In lattice NRQCD, the heavy quark is described in 
terms of a two-component (in spin space) field i/j(x). The 
Dirac field h(x) is related to ip(x) (and the antiquark field 
?/>(#)) by a unitary Foldy-Wouthuysen transformation [0, 



(5) 



This transformation decouples the upper and lower com- 
ponents of the Dirac field, thereby separating the quark 
field from the antiquark field. To facilitate the identifi- 
cation of lattice NRQCD operators capable of matching 
Eq. (^), we similarly transform the external state spinor 
Uh(p) into a nonrclativistic Pauli spinor: 



Uh{p) = 



where 



)(P) 



u Q (p)+0(l/M 2 ), (6) 



V 



and Uq is a two-component external state spinor depend- 
ing only on the spin of the heavy quark. Note that we are 
working in the Dirac-Pauli representation. Using Eq. (|^), 
the relation 70 uq (p) — uq (p) , and the Dirac equation 
for the light quark uJp')p'o = — Wg(p')"y-p'7o, the spatial 
components of Eq. (0) may be written 



[q(p')\J k \h(p)) 



QCD 



(1) 



-0(a 2 ,l/M 2 ) 



(4) 
k 



with 



««(?') r fc U Q(P)> 

-iUq( 

iu q (p' 



iUqip') r fe^J u Q(p)i 



-^-70 r fe u Q (p), 
-™ q (p') & r J u Q ( P ). 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
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The coefficients in Eq. (||) are given by 

Vo = ai = 1 + a s B , 
Vi = Vo, 

T] 2 = 2a 3 = a s B 2 , 
T] 3 = 2a 2 = a s B 3 , 
r/4 = 2 04 = a s B4, 



(14) 



where 



B a = — 
B 2 
B 3 
Bi 



3i 4-t 
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Having obtained the 1 /M expansion of the above scat- 
tering amplitude in continuum QCD, the next step is to 
identify operators in the lattice effective theory that can 
reproduce the terms in this expansion. An inspection of 
Eq. (||) suggests immediately that matrix elements of the 
following five lattice operators should be considered: 



^ifc,lat( a; ) 


= q{x) r fe 7 Q(x), 


(16) 


J (1) (x) 




(17) 


/ (2) (x) 


= ^eW 7-^70 r{Q( x ), 


(18) 




= ^q{x)Y J Q V k Q{x) 


(19) 


J (4) (x) 




(20) 



where q(x) is now the light quark field in the lattice the- 
ory, Mq is the bare heavy quark mass, and Q{x) is related 
to the heavy quark field ip(x) in lattice NRQCD by 



Q(x) 



ijj(x) 




(21) 



We use the bare quark mass Mg in the above definitions 
since it is the natural choice to use in lattice simulations 
and because the pole mass M is not well defined beyond 
perturbation theory. Definitions of the lattice derivatives 
are given in the next section. 

For the temporal components of the currents Jo = Vo 
or Aq, the leading 1/M behavior is given by 



[q(p')\Jo\h(p)) 



QCD 



>(2) 



+0(a 2 s ,l/M 2 ), (22) 



n^ = u q ( P ') r J u Q (p), 



(i) 



(2) 



-iu q (p') rQ^j U Q (p), 

iu q (p') "|^-7o Tq 7 u Q {p). 



n, 



The coefficients in Eq. (E2J) are given by 
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Vq = («1 + a 2) 



l + a s Bq, 



(a-i - a 2 ) = 1 + a s B{, 
r\\ = 2 (a 3 + a 4 + a 5 ) = a s B\, 



(26) 



where 
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Again, it is straightforward to identify operators in the 
lattice effective theory which can reproduce these terms: 



^0,lat( 2: ) 
"A), lit ( X ) 



1 



2M 
-1 

2Mq 



q(x) r^.VQ(i), 



(x) 7 -^7or o 7 Q(x). 



(28) 
(29) 

(30) 



These operators were previously used in the expansion of 
the temporal component of the axial-vector heavy-light 
current in Refs. 



III. MIXING MATRIX 

Having identified the necessary operators in the lattice- 
regulated effective field theory, we next need to calculate 
the mixing matrix of these operators defined by 



( i) 

where fljj [' are the five fl k 



with 



W(a 2 , 1/M 2 a 2 a s a/M), (31) 

{J) defined in Eqs. (|)-(|l3|) for 

/i = k = 1,2,3 and the three fig defined in Eqs. (^~~ 

7 j 



( |25| ) for /j = 0. The determination of the Z^ factors 
is the most difficult aspect of the current construction. 
These factors are determined numerically using lattice 
perturbation theory. The same methods used in Ref. Q 
are applied here. We also use the same lattice action 
as in Ref. H, but for the convenience of the reader, we 
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restate these actions below. The Feynman rules for our 
heavy and light quark lattice action and further details 
of lattice perturbation theory are given in Refs. HQ- 

For the heavy quark, the following NRQCD action den- 
sity on the lattice Q is used: 

clCnrqcd = ip^(x) ip(x) 



-il)\x+at)(l- 



,n)Il \f aH a \ n u\{x) 



x 1 



aH Cl 
2n 



2 J \ 2n 
aSH 



where 



H 



A( 2 ) 
2Mo' 



SH = -c B 



2Mq 



B, 



Uq 

i>(x), (32) 

(33) 
(34) 



U^x) are standard gauge-field link variables, n is an in- 
teger, <jj are the Pauli matrices, uq is the mean link pa- 
rameter introduced by the tadpole improvement proce- 
dure flCj ], and the QCD coupling g is related to a s by 
a s ={f /(AiT). At tree level, eg = 1; the one-loop con- 
tribution to cb is an 0(a 2 ) effect in our mixing matrix 
calculation and hence can be ignored here. The chromo- 
magnetic field is given by Bj{x) — — \(.ji m Fi m {x) 1 where 
the Hermitian and tracclcss field strength tensor F^(x) 
is defined at the sites of the lattice in terms of clover-leaf 
operators: 



■F i u,v(x) 



2a 2 g 



tonv(x) = —j ^ U a (x)Up(x+a a ) 

x U- a (x + a a +af3)U-p(x + a{3), 



(35) 



with {(a, 0)}^ = {(/i, v), (v, -p,), (-^, -u), {-v, fi)} for 
/.t v. The lattice derivatives in our action and in the 
current operators are given by 



(36) 
(37) 



a\/^0(x) = —[U^xpix+a^) 

-Ufa-a^O^-a^)], 

-O^x-a^U^x-a^)], 

3 

i J A (J) 0(i) = ^2fa l [U k (x)0(x+a k ) 
fe=i 

+Ul(x-a k )0(x-a k )} - 20(x)), (38) 

3 



+Ul(x-a lt )0(x-a lt )] - 20(x)), (39) 

where 0{x) is an operator defined at lattice site x with 
appropriate color structure. 

For the light quarks, we use the clover action 0], 



aCught = qfq- a-qV^q + m qq 



i9a J^2 Q^nvF/iu q, 



(40) 



where J/ = 5^7/jV^, m o i s the bare light quark mass, 
v^v = \ [7/i, 7«]j and we set the Wilson parameter r = 1. 
The one-loop correction to the clover coefficient is an 
O(a^) effect in our matching calculation and is neglected 
here. 

At one-loop order, the mixing matrix elements may be 
written Ol 



{[h(c q 



+C Q )+C m (l-S i0 ) 



}■ 



(41) 



where C q and Cq are the contributions from the light - 
and heavy-quark external leg corrections (that is, from 
wave function renormalization factors), and denote 
the contributions from the vertex corrections. Our use of 
an on-shell renormalization scheme with lattice operators 
defined in terms of the bare mass Mq is responsible for 
the term proportional to C m , where 

M = [l + a s C m ] M + 0(a 2 ). (42) 

Note that although the current operators Jy] at are de- 
fined using Mo, the pole mass M appears in Q^. Within 
the context of finite-order perturbation theory, the pole 
mass is an observable quantity. It is gauge invariant and 
rcgularization scheme independent. Hence, it is the nat- 
ural heavy-quark mass definition to use when matching 
amplitudes in different schemes. However, in reality, a 
confined quark has no pole mass. Once the matching 
coefficients Cj are obtained in terms of aM, it is then 
necessary to eliminate any dependence on M in favor of 
some quark mass which is defined beyond perturbation 
theory. Clearly, the mass parameter Mq appearing in the 
lattice action is the most suitable and convenient choice. 
The relationship between M and Mo given in Eq. 
is then used to re-express the Cj 7 coefficients in terms of 

Mq. 

The factors in Eq. (HJ) may be further decomposed: 
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37 



C q = C q + ^lnaX + Cj 1 , 



C Q = C q - — InaX, 



Cm — Cn 



a 



TI 



C J - =C J - +( 



J,TI 



(43) 
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where C q , Cq, C m , and (fj are infrared finite and inde- 

J IR, 

pendent of the tadpole improvement factor uq, and Qj 

J TI 

and Qj contain the infrared divergences and tadpole 
improvement contributions, respectively, from the vertex 
corrections. Contributions to C q and C m from the tad- 
pole improvement counterterms are denoted by C q l and 
C^ 1 , respectively. The tadpole improvement terms are 



C, 



TI 



.(2) 



r TI - -7/ 



(2) 




1 - 



-J,TI 



,( 2 ); 



2 n aMo 



(44) 



where uq = 1 



,( 2 ) 



O(a^). For the usual plaquette 



definition u = (^TWo) 1 / 4 in the Wilson gluonic action, 



,(2) 



7r/3. For massless clover quarks, C, 



Values for Cq and C m are given in Ref. 



1.030. 



IV. MATCHING 



To complete the operator construction for J/., we in- 
vert Eq. © to get n<?> = E i (Z J )j l 1 (jjgjiat, substi- 
tute this into Eq. (||), peel off the dependence on the 
external states, and use Eqs. (|l4|), (15), (fli]), (42), and 
([43|) to obtain 



J fe = (1 + a s 
\ + a s 



Bn — Cn n — 
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-|ln(aM ) + i 



and 



C Qq = 9(^9 + C'J 1 + Cq)' 



CQqm — Cq<3 ~\~ Cm + C, 



TI 



(46) 

(47) 
(48) 



As expected, the infrared divergences from the various 
terms in the expansion coefficients cancel. Results for 
rj k and r^ k for various values of aMg are listed in Ta- 
bles I and [II], respectively. 

For the temporal component of the vector current, we 
substitute the inverted Eq. ([|l]) into Eq. (^||), remove 
the dependence on the external states, and use Eqs. (|2£ 
@, (Hh, @, and (ph to obtain 
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Again, all infrared divergences cancel. Results for for 



various aMo values are given in Table III 
and 



From Eqs. ([45j) and (|49j), one easily obtains the match- 
ing coefficients Cj of Eq. ([!]) which expresses the heavy- 
light currents in terms of operators suitable for numerical 
simulations using lattice NRQCD. These coefficients can 
be written in the form 



Ci(a s ,aM ) 



l + a s pf(aM ) + O(a 2 ), (i = 0,1), 
a s p{ (aM ) + 0(a 2 ), (i > 2). 

(51) 



The values for p^ h , pf k , and p\° for various values of 
aMo are listed in Tables |y|, and VI, respectively. 
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V. 0(aa s ) CORRECTIONS 

In this section, discretization corrections in the static 
limit aMo — > 00 are discussed. To study such correc- 
tions, we calculated all rf for several large values of the 
heavy quark mass. With the exception of r/, all of the 
t/ tended tamely to finite values in the static limit. How- 
ever, as aMo became large, we found that the magnitude 
of t/ began to grow linearly with aMo: 



2 aM ^ 1 - 2aM C'Lc + constant, 



(52) 
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where the disc are constants independent of aMo- We 
determined the disc computing t/ for aMo — 
10.0,25.0,100.0,400.0,1000.0 and 5000.0. We then fit 
the results for — 7"//(2aMo) to a quadratic polynomial 
in l/(aMo) and obtained the following limiting values: 



CZc= 1-029(2), (53) 

C& c = L031 (l), (54) 

Ctc = -0.063(1), (55) 

CaL = -0-063(1). (56) 



The discretization factor £^° c for the temporal compo- 
nent of the axial- vector current was previously calculated 
in Ref. [0] , but its determination was not done very accu- 
rately. We have recalculated Cdisc here to a much higher- 
precision. 

Since the contributions proportional to disc are P urer y 
0(aa s ) discretization corrections, one may absorb these 
terms into the lattice current operators. Improved cur- 
rent operators can be defined using 

■>SLt(s) = J St(*) + CLc ^ffW. (57) 
^lt(^) = ^LW 5 (i>0), (58) 
where 

= -aq(x) 1 .^ 1 vVlQ{x). (59) 
The mixing of the improved operators is given by 

(9(p')I^Ltl%)}lat = E 

3 

+0{a 2 s ,l/M 2 ,a 2 ,a s a/M), (60) 

where 

Z£ = Zi + 2a s aM CLc hi , (61) 
£/i = ^, (i>0). (62) 

Then the coefficients of these operators are 

Ci = Ci-2aM disc**, (63) 
Cf = Cf, (i^2), (64) 

and the heavy-light currents are given by 

^ = E^/ J Sat + °( a l a " l l M l <*sa/M Q ). (65) 

i 

Terms which grow linearly with aMo are no longer 
present in C'j and Z(=. Note that in the lattice NRQCD 
approach, discretization and relativistic corrections are 
intertwined since 0(a) and 0(1/ M) interactions are 
treated as equally important. 



VI. LARGE LOGARITHMS 



Another feature of the Cj matching coefficients is the 
presence of ln(aMo) terms. The simulations of heavy- 
light systems carried out to date have used bare heavy 
quark mass values in the range aMo ~ 1-6 — 4.0 where 
these logarithms are not large. However, in the large 
aMo limit, these logarithms must be treated with care. 
The logarithms appearing in the coefficients for j > 
are tamed by the l/(aMo) factors in their corresponding 
current operators. However, the logarithm appearing in 
the Cq coefficient becomes problematical and must be 
dealt with using the renormalization group. 

In matching calculations between QCD and various 
continuum effective theories, one usually encounters sim- 
ilar logarithms of the form \n(M/fjb), where /x is some 
scale introduced by the renormalization procedure. Such 
logarithms are summed using the renormalization group 
(RG) equations which follow from the requirement that 
physical quantities must not depend on fi. Since \x ap- 
pears only inside the logarithms, one ends up with simple 
anomalous-dimension matrices and RG equations which 
can be solved straightforwardly. The situation is more 
complicated in the present calculation. The role of /x is 
taken over by the inverse lattice spacing 1/a, and a ap- 
pears not only in the logarithms, but also in other places, 
such as the Tj 7 which are complicated functions of aM . 
Furthermore, the ultraviolet cutoff 1/a is an integral part 
of our effective theory and we cannot take a — > 0. 

As discussed in Ref. Q, the observation that the left- 
hand side of Eq. ([!]) is independent of the lattice spacing 
can be exploited to derive an RG equation for the Cj 
coefficients. This equation describes the change in the 
Cj as the lattice spacing is varied. Collecting the co- 
efficients into a vector, the RG equation may be written 

(a^- + ( 7 J ) tr ^ S-J = 0, (66) 

where the anomalous dimension matrix is given by 

^(a s ,aM Q ) = ]T U-^zA (Z%\ (67) 

k ^ ' 

The right-hand side of Eq. (|6^) is a complicated function 
of aMo- Once 7^ is determined numerically for a large 
range of aMo values, Eq. ( |66| ) can be solved by numerical 
methods. Our primary concern was the determination of 
the matching coefficients for values of aMo relevant 
for simulations of B and D mesons, and for such val- 
ues, RG improvement was not needed. Hence, we have 
not attempted to obtain the entire anomalous dimension 
matrices for our current operators. 
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VII. SUMMARY 

In summary, the spatial components of the heavy- light 
axial- vector current Ak and all components of the heavy- 
light vector current were expressed in terms of lat- 
tice operators suitable for use in simulations of B and 
D mesons. In the lattice theory, the heavy quarks were 
treated using the NRQCD formulation, the light quarks 
were described by the tadpole-improved clover action, 
and the standard Wilson action was used for the gluons. 
The light quarks were treated as massless. The expan- 
sions were carried out to 0(1/ M) by matching appropri- 
ate scattering amplitudes to one-loop order in perturba- 
tion theory. We found 

J»=J2 C I J Sat + 0(a 2 s ,a 2 , l/Ml a s a/M ), (68) 

3=0 



where N^ h = Ny k = 5 and Ny 
operators are given by 



3. The lattice current 



(0) (x 



J (1) (r 



t(2) 



(x 



J (4) (x 



q(x)TiQ(x) 



q(x)r' 7 .VQ(x) 



1 q(x) ltfV fc Q(x) 



2M 
1 

2Mq 



(69) 
(70) 

(71) 

(72) 

(73) 



with T% = 7 „ and r£ 



ci 



757/j: an< i the coefficients are 
1 + a s P j{aM Q ) + 0(a 2 s ), (i = 0, 1), 



a s pi(aM ) + O(a 2 s ), (i > 2), 



(74) 



where values of the pf (aM ) are listed in Tables IV VI. 
The currents can also be expressed in terms of improved 
current operators J^] at (x) as shown in Eq. (|65|). 

This completes the matching calculation through 
0(a s /M) and 0(aa s ) for all components of the vector 
and axial-vector heavy-light currents. Our matching co- 
efficients have already been applied in leptonic B and B* 
meson decays to extract the fps and fy decay constants 
p2|-p^[. They are also relevant for studies of B — * n or 
p semileptonic decays. In this article, we presented re- 
sults only for the simple NRQCD action of Eqs. (|3l)-((34|)- 



Matching coefficients for other NRQCD actions which 
have appeared in the literature and for different values of 
(aMo, n) are also available. For example, Ref. Jl^] used a 
slightly different action, and an action with higher-order 
improvement terms was employed in Refs. [ [U3||lE| ]. In all 
cases, we find that in the useful range 1.0 < aMo < 10.0, 
the one-loop coefficients exhibit only a mild dependence 
on aM and do not become particularly large. 
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TABLE I. Values of the coefficients r i h corresponding to the spatial components of 
the vector current Vk for various values of the bare heavy-quark mass aMo and NRQCD 
stability parameter n. Uncertainties in the determinations of these parameters due to the 
use of Monte Carlo integration are included. 



aM 


n 


v k 
T o 


Vk 


T 2 


Vk 
T 3 


Vk 
T 4 


fo.o 


1 


0.8189(1) 


— 0.588(3) 


-13.93(1) 


-0.566(4) 


-0.940(5 


7.0 


1 


0.7880(1) 


— 0.640(2) 


-8.150(7) 


-0.508(3) 


-0.909(3 


4.0 


1 


0.7221(1) 


-0.757(5) 


-2.809(7) 


-0.385(7) 


-0.836(7 


4.0 


2 


0.7340(1) 


-0.742(4) 


-3.084(8) 


-0.422(6) 


-0.842(7 


3.5 


2 


0.7176(2) 


-0.775(4) 


-2.310(5) 


-0.392(6) 


-0.817(6 


3.0 


2 


0.6982(1) 


-0.821(4) 


-1.594(4) 


-0.351(5) 


-0.791(5 


2.7 


2 


0.6840(2) 


-0.856(3) 


-1.196(4) 


-0.323(5) 


-0.767(5 


2.5 


2 


0.6732(2) 


-0.877(3) 


-0.955(4) 


-0.309(5) 


-0.750(4 


2.0 


2 


0.6410(1) 


-0.962(3) 


-0.427(3) 


-0.254(4) 


-0.693(3 


1.7 


2 


0.6170(2) 


-1.027(2) 


-0.188(2) 


-0.218(3) 


-0.647(3 


1.6 


2 


0.6075(2) 


-1.051(2) 


-0.123(2) 


-0.205(3) 


-0.627(2 


1.2 


3 


0.5768(2) 


-1.175(2) 


-0.081(2) 


-0.150(3) 


-0.534(2 


1.0 


4 


0.5559(2) 


-1.265(2) 


-0.148(1) 


-0.108(3) 


-0.464(2 


0.8 


5 


0.5249(2) 


-1.387(2) 


-0.288(1) 


-0.049(3) 


-0.363(1 



TABLE II. Values of the coefficients T i * corresponding to the spatial components of 
the axial- vector current Ah, similar to Table M. 



aM 


n 


Ak 
' 


Ah 




T 3 


T Ak 


10.0 


1 


0.2634(1) 


-1.098(4) 


2.624(18) 


0.611(6) 


-0.7291(2 


7.0 


1 


0.2685(1) 


-1.145(3) 


2.137(12) 


0.636(4) 


-0.7017(2 


4.0 


1 


0.2817(1) 


-1.254(2) 


1.494(6) 


0.691(2) 


-0.6401(2 


4.0 


2 


0.2742(1) 


-1.197(2) 


1.475(6) 


0.590(2) 


-0.6407(2 


3.5 


2 


0.2766(1) 


-1.225(2) 


1.339(5) 


0.595(2) 


-0.6219(2 


3.0 


2 


0.2793(1) 


-1.259(1) 


1.168(4) 


0.598(2) 


-0.5977(2 


2.7 


2 


0.2812(1) 


-1.284(1) 


1.065(3) 


0.598(2) 


-0.5796(2 


2.5 


2 


0.2824(1) 


-1.303(1) 


0.985(3) 


0.600(2) 


-0.5654(2 


2.0 


2 


0.2853(1) 


-1.365(1) 


0.750(2) 


0.599(2) 


-0.5196(2 


1.7 


2 


0.2862(1) 


-1.416(1) 


0.582(2) 


0.595(2) 


-0.4816(2 


1.6 


2 


0.2860(1) 


-1.436(1) 


0.513(2) 


0.591(2) 


-0.4661(2 


1.2 


3 


0.2725(1) 


-1.510(2) 


0.188(1) 


0.537(2) 


-0.3823(2 


1.0 


4 


0.2598(1) 


-1.567(2) 


-0.025(1) 


0.511(2) 


-0.3190(2 


0.8 


5 


0.2407(1) 


-1.658(2) 


-0.291(1) 


0.500(2) 


-0.2290(2 
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TABLE III. Values of the coefficients r 4 correspond- 
ing to the temporal component of the vector current Vo, 
similar to Table fl. 



aM 


n 


r 




T 2 


10.0 


1 


0.3549(2) 


-0.541(8) 


1.706(19) 


7.0 


1 


0.3961(2) 


-0.572(6) 


1.339(12) 


4.0 


1 


0.4954(2) 


-0.643(3) 


0.982(12) 


4.0 


2 


0.4868(1) 


-0.695(3) 


0.969(6) 


3.5 


2 


0.5161(1) 


-0.718(3) 


0.907(5) 


3.0 


2 


0.5542(1) 


-0.752(3) 


0.844(4) 


2.7 


2 


0.5832(2) 


-0.774(2) 


0.815(4) 


2.5 


2 


0.6058(1) 


-0.794(2) 


0.792(3) 


2.0 


2 


0.6795(2) 


-0.853(2) 


0.744(3) 


1.7 


2 


0.7417(2) 


-0.899(2) 


0.718(2) 


1.6 


2 


0.7666(2) 


-0.920(2) 


0.707(2) 


1.2 


3 


0.8832(2) 


-1.034(1) 


0.684(1) 


1.0 


4 


0.9674(2) 


-1.109(1) 


0.682(1) 


0.8 


5 


1.0853(2) 


-1.202(1) 


0.696(1) 



TABLE IV. Values of the coefficients p- k defined in Eq. (|3l|) corresponding to the 
spatial components of the vector current for various values of the bare heavy-quark mass 
aMo and NRQCD stability parameter n. Uncertainties in the determinations of these 
parameters due to the use of Monte Carlo integration are included. 



aM 


n 


v k 
Po 


v k 
Pi 


P 2 


v k 
Pa 


P4 


10.0 


1 


-0.5051(3) 


-0.128(5) 


16.968(11) 


0.991(4) 


0.069(5 


7.0 


1 


-0.5441(2) 


-0.202(4) 


10.733(7) 


0.932(3) 


0.190(3 


4.0 


1 


-0.5496(2) 


-0.346(5) 


4.680(7) 


0.810(7) 


0.354(7 


4.0 


2 


-0.5744(2) 


-0.382(5) 


4.955(8) 


0.846(6) 


0.360(7 


3.5 


2 


-0.5694(2) 


-0.421(5) 


4.011(5) 


0.816(6) 


0.392(6 


3.0 


2 


-0.5585(2) 


-0.443(5) 


3.099(4) 


0.775(5) 


0.431(5 


2.7 


2 


-0.5472(2) 


-0.470(4) 


2.566(4) 


0.747(5) 


0.451(5 


2.5 


2 


-0.5366(2) 


-0.495(4) 


2.228(4) 


0.734(5) 


0.467(4 


2.0 


2 


-0.4958(2) 


-0.521(3) 


1.415(3) 


0.678(4) 


0.505(3 


1.7 


2 


-0.4561(3) 


-0.551(3) 


0.970(2) 


0.642(3) 


0.528(3 


1.6 


2 


-0.4391(3) 


-0.564(3) 


0.828(2) 


0.630(3) 


0.534(2 


1.2 


3 


-0.3679(3) 


-0.609(3) 


0.419(2) 


0.574(3) 


0.563(2 


1.0 


4 


-0.3018(4) 


-0.606(3) 


0.254(1) 


0.532(3) 


0.570(2 


0.8 


5 


-0.1818(5) 


-0.597(4) 


0.110(1) 


0.473(3) 


0.564(1 
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TABLE V. Values of the coefficients pf k corresponding to the spatial components of 
the axial- vector current, similar to Table IV. 



aM 


n 


Po 


A h 
Pi 


At 
P2 


Ai 
P 3 


Pa 


10.0 


1 


0.0504(2) 


0.383(6) 


0.413(18) 


-0.186(6) 


-0.1420(2 


7.0 


1 


-0.0246(2) 


0.303(4) 


0.447(12) 


-0.211(4) 


-0.0180(2 


4.0 


1 


-0.1093(2) 


0.151(3) 


0.377(6) 


-0.267(2) 


0.1578(2 


4.0 


2 


-0.1146(2) 


0.073(3) 


0.397(6) 


-0.166(2) 


0.1584(2 


3.5 


2 


-0.1285(1) 


0.028(3) 


0.362(5) 


-0.171(2) 


0.1963(2 


3.0 


2 


-0.1396(1) 


-0.005(3) 


0.337(4) 


-0.174(2) 


0.2375(2 


2.7 


2 


-0.1444(2) 


-0.043(3) 


0.305(3) 


-0.174(2) 


0.2641(2 


2.5 


2 


-0.1458(2) 


-0.070(3) 


0.288(3) 


-0.176(2) 


0.2826(2 


2.0 


2 


-0.1401(2) 


-0.117(3) 


0.238(2) 


-0.174(2) 


0.3316(2 


1.7 


2 


-0.1253(2) 


-0.161(3) 


0.200(2) 


-0.171(2) 


0.3625(2 


1.6 


2 


-0.1177(2) 


-0.180(3) 


0.191(2) 


-0.166(2) 


0.3727(2 


1.2 


3 


-0.0636(3) 


-0.275(3) 


0.151(1) 


-0.112(2) 


0.4110(2 


1.0 


4 


-0.0058(3) 


-0.304(3) 


0.131(1) 


-0.087(2) 


0.4251(2 


0.8 


5 


0.1023(5) 


-0.326(4) 


0.113(1) 


-0.076(2) 


0.4298(2 



TABLE VI. Values of the coefficients p-° correspond- 
ing to the temporal component of the vector current, sim- 
ilar to Table |jv]. 



aM 


n 


Vh 

Po 


Vh 
Pi 


Vh 
Pi 


10.0 


1 


0.1712(3) 


-0.387(9) 


-0.433(19) 


7.0 


1 


0.0599(2) 


-0.482(7) 


-0.065(12) 


4.0 


1 


-0.1107(3) 


-0.672(4) 


0.291(12) 


4.0 


2 


-0.1150(2) 


-0.642(5) 


0.305(6) 


3.5 


2 


-0.1557(2) 


-0.691(4) 


0.367(5) 


3.0 


2 


-0.2023(2) 


-0.725(4) 


0.430(4) 


2.7 


2 


-0.2341(2) 


-0.765(3) 


0.458(4) 


2.5 


2 


-0.2570(2) 


-0.791(3) 


0.481(3) 


2.0 


2 


-0.3221(2) 


-0.842(3) 


0.529(3) 


1.7 


2 


-0.3686(3) 


-0.891(3) 


0.555(2) 


1.6 


2 


-0.3861(3) 


-0.907(3) 


0.566(2) 


1.2 


3 


-0.4621(3) 


-0.963(3) 


0.589(1) 


1.0 


4 


-0.5012(4) 


-0.974(3) 


0.591(1) 


0.8 


5 


-0.5300(5) 


-0.994(3) 


0.577(1) 
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